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ABSTRACT 

We study the corrections to the determination of sin(2/5) from the time de- 
pendent CP asymmetry of J /'^ Kg which arise in the standard model. 
Although a precise prediction of these corrections is not possible we find 
that they are indeed extremely small, of the order of less than a per mil of 
the observed value. This means in turn that any deviation visible at the B 
factories will be a clear signal for new physics. 
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1 Introduction 



The measurement of the mixing- induced CP asymmetry in the so-called gold-plated mode 
J/'^ Ks is becoming a precision measurement jT]. Currently the relative uncertainty is 
at the level of five percent and is going to decrease further in the next few years as more data 
from the B factories are analyzed. 

From the theoretical side the time dependent CP asymmetry in this channel is related 
in the standard model to sin(2/3) in a very clean way i.e. it is not plagued by hadronic 
uncertainties, at least at the level of current experimental precision. However, the precision of 
the experimental data will increase further, making this CP asymmetry an interesting probe 
for new physics. 

Possible new physics effects have been discussed in some detail in a generic framework in j3] 
where both the charged and the neutral modes oi B ^ J /^! K have to be taken into account. In 
|1] certain observables have been defined which are sensitive to different aspects of new physics. 

However, in order to quantitatively pin down a possible new physics effect and to assess 
the reach to new physics, the small standard-model contributions have to be under control. 
Although estimates of these effects have been given some time ago ^5J, it is worthwhile to 
reconsider these estimates motivated by the experimental precision to be expected soon. 

In this paper we shall try to get a quantitative estimate for the time-dependent CP asym- 
metry in this channel beyond the simple relations (in the convention used in below) 

Cb^j/^Ks = Sb^j/^ Ks = sin(2/3) . 

The corrections to these relations originate from two sources. The first source is from the 
corrections to the /S.B = ±2 part of the Hamiltonian. In the present paper we compute these 
contributions systematically in an effective field theory by subsequently integrating out heavy 
degrees of freedom down to a small hadronic scale. It is worthwhile to note that many elements 
of our calculation are similar to the calculation of the lifetime difference in the system of neutral 
B mesons which are performed in an effective field theory framework in [0]. 

The second source is the decay amplitude itself which has small contributions carrying a 
different CP phase compared to the leading piece. The latter is much harder to estimate and 
we shall refer to well known methods. 

A full discussion of the standard-model effects would also require to include the effects from 
CP violation in the Kaon system. We shall not discuss these effects in the present paper, we 
focus completely on the effects coming from the B system 

In the next sections we set up the calculation of these two contributions and discuss the 
methods of calculation. Finally we summarize our results and conclude. 

2 Basic relations 

CP asymmetries are measured at the B factories by detecting the decay products of the coherent 
B^B^ pair. One of them is identified as a flavor state using e.g. a leptonic tagging mode while 
the decay of the other into a CP eigenstate is observed. Usually, in the calculation of these 
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rates a possible CP violation stemming from the mixing on the tagging side or in the tagging 
decay is neglected which results in the well-known formula for the decay J/'^Kg 



a. 



B^JI-i! Ks 
CP 



it) 



Cb~.j/^Ks cos(Amt) - SB^j/'jfKs sin(Amt) 
cosh(Ar t/2) + De^j/^f Ks sinh( AF t/2) 



with 
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2 Re [A] 

1 + |A|2 
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(3) 



Here p and q are the parameters characterizing the mixing in the K- and -B-system, while 
A{B^ — i> J I'^Ks) is the amplitude for the decay into the CP eigenstate J /^Ks- 

Our aim is to investigate all effects correcting the leading contribution and therefore we 
take also into account the influence of the mixing in the B system on the tagging amplitude 
assuming that there is no CP violation in the tagging decay (cf. f?). Defining 



l^tagl 



the CP asymmetry becomes 
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B^J/^Ks 
CP 



it) + - 
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sin (Amt) sm\2P). 



(4) 



(5) 



These relations above are true in general and take into account a possible width difference, 
a possible direct CP asymmetry as well as the mixing effects on the tagging. The lifetime 
difference AF as well as 1 — |Ap and e are small and the leading term is obtained by neglecting 
these quantities. Furthermore, since the weak phase of the decay amplitude of the leading 
contribution vanishes in the standard convention, S measures the weak phase of the AB = 2 
contribution to the effective Hamiltonian which is - again to leading order - simply 2f3. 

We make use of the general relation 



B^J/^ Ks 



B-^J/'S'Ks 



B^J/'S'Ks 



)'■ 



(6) 



which allows us to replace Db^j/^Ks by its leading order expression Db^j/^Ks = cos(2/?), 
since it is multiplied by the small quantity ATt in the expansion. Note that typically t is of 
the order of the B meson lifetime r and we have AF r <^ 1 . 



3 Corrections to the mixing 

The phenomenon of mixing between B and B is due to the box diagrams with a double W 
exchange in the full electroweak theory. These diagrams have been evaluated some time ago 
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in the full standard model including the quark masses of all quarks in the loop . After GIM 
cancellation, the leading term is due to the top quark, giving rise to a contribution of the 
order {mt/My/Y' with the weak phase 2(3. Subleading terms are either of the order (mt/Miy)^ 
which again carry the phase 2(3 and thus will not contribute to a modification of S", or of order 

/Mwf which carry a different weak phase and hence will yield a correction to S. 

Instead of calculating the box diagrams in the full theory one may also use an effective 
theory picture which we will do in this paper. The approach is similar to the one used for 
D — D mixing jS^QflU j. and also for the AS" = 2 Hamiltonian in fTD . The advantage is that one 
may express the contributions to the mixing phase in terms of matrix elements of certain local 
operators, which can be estimated in factorization. Furthermore, one can resum systematically 
large logarithms using the renormalization group. 

The first step is the usual one, in which the W boson and the top quark are integrated out 
at a common scale ^ rrit ^ My/. The result is the well known Ai? = 2 operator of dimension 
six, which is 

H^E=^ = ^>^lmlC^{xt)Qo with Q, = {hLlMi^LYdL) (7) 
where \q = V*f^Vqd, Xt = mf/M^ and 

4 - llxt + x1 ?,xl\Yixt , . 

^{\-XtY 2(1 -a;i)'^ 



is obtained from the usual Inami-Lim function |T2] . Note that we extracted explicitly a factor 
vn\ which is a reminiscent of the GIM mechanism. 

However, (0) is not the only contribution to Ai? = 2. Other contributions originate from 
two insertions of Ai? = 1 operators 

T^^- = -1 1 d^xT {E^i-\x)E^i-\^)\ , (9) 

where the relevant operators are 

AC 

H^^=^ = [V:,VUhl,CL) (cLl'dL) + V:,VUbLl,CL) iuLl^dL) (10) 

+V*^ycd{hLl^,UL){cLl^dL) + V*^Vud{hLl^,UL){uLl^dL)\ 

which yield non-local contributions at the scale /i ~ Mw- 

Furthermore, the matching at the scale Mw yields - aside from the above dimension-6 op- 
erator appearing in (0) - dimension-8 contributions involving two covariant derivatives. These 
have been partially calculated in ^3]; however, we do not need these contributions for our 
purposes, since these pieces are again dominated by the top quark and thus have the same 
weak phase as the leading part. 

In addition, keeping a non-zero charm mass, the matching calculation yields an operator 
of the form m^(6i7^(ii)(6i,7^(ij^) which we shall treat as a dimension-8 operator as well. The 
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matching of this operator involves the calculation of the box diagrams keeping the charm mass 
non-zero. 

Lowering the scale turns the high momentum part of the non-local contributios into local 
operators. This effect is described by a renormalization-group mixing of the non-local operators 
into local ones. We define the non-local operators as 

Ti = ~J d''xT[{h^^CL){x){cLrdL){x){hrcLm^^^^ (11a) 



T2 = -7: d^xT[{bL'y^,CL){x){uLrdL){x){bLYuL)mc^^^ 



d^X T [{h-f^UL) {X) {uLYdL) {X) {bLYuL){0) {uLl.dL) (0)] 



(lib) 

(lie) 



and find that these operators mix already at order (0:5)° into local AS = 2 operators of 
dimension 8^ 



= U{hL^^dL){hLl^dL) 

Q2 = d''d''{bL7MibL7.dL) 
Qs - mlipLlM^Ll^dL). 

This order (q;s)° mixing happens via the diagrams 

& u.c d 




d ' h 

The anomalous dimension matrix at order (q;s)° is 
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(12) 



(13) 



and the operators are gathered into O — (Qi, (52, Qs, ^i, ^2, T^Y . 

^If the complete matching at /x = Mw would be performed, many more dim-8 operators would be induced; 

however, their coefficients would not be enhanced by a large logarithm; large logarithms at order (o-sY' are only 
induced by the mixing with the non-local terms, which induces only Qi and Q^- 
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The initial conditions for the renormahzation-group running are from matching. As dis- 
cussed above, we do not exphcitly need the matching conditions for Qi and Q2, so we get 



C^iMw,Xt) = -^Ci2(xi), C2iMw,Xt) = -^C[,ixt), 

Cs{Mw,xt) = ^ {2X,XtC3{xt) + Xl) , (14) 



C^iMw) = Xl C,{Mw)=2X,X^, Ce{Mw) = X 



where we have used the tree-level matching of the Ai? = 1 operators for the non-local terms. 
The function 

C3(.,)^.n.,-^-i£ill|L (15) 

is derived from the Inami-Lim function; C12 and €[2 are functions of Xt which could be obtained 
from the matching of the dim-8 operators with two derivatives. We do not give these functions 
here, since they are not needed in the following. 

The next step is to perform the renormalization-group running using the renormalization- 
group functions at order (0;^)° 

^^^^-l''^c{^i) = Q (16) 

where C = {Ci, C2, C3, C^, Cq) are the coefficients of the operators. 

Running down to the scale of the b quark mass we switch at /i ~ rrib again to another 
effective field theory in which the b quark becomes static. At this scale one has to replace the 
derivatives acting on the b quark field by 



idf,b — > {rribV^ + idfj)hy (17) 

where is the static b quark field moving with velocity v. Keeping only the leading term in 
the l/rrif, expansion we have to match onto the operators 

P3= rnlChi'-llMChi-lrdL), 



where h[^l~^ denotes the static quark/antiquark field, which have become completely different 
fields in the static limit. 

Performing the matching and the renormalization-group running we may evolve down to 
scales around the charm-quark mass. At such a low scale all contributions become local, once 
the up-quark mass is neglected. In fact the dim-8 operators of the generic structure 
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with D being a derivative inserted in any possible way, will yield contributions proportional to 
^QCD- Compared to the pieces considered here, these are suppressed by a factor K^^Qj^/ml and 
thus may be neglected. 

Inserting back all the CKM factors, we make use of CKM unitarity, i.e. the fact that 
Am + Ac + At = 0. This separates the running from CKM factors, which means that each weak 
phase is multiplied by a renormalization-group invariant. Thus the A5 = 2 contributions at 
the scale /i ~ can be expressed in terms of the local effective Hamiltonian 



AB= 
eff 



4^' 



2A? 



A? 



C'uixt) - In 



Ci2{xt) - In 



A + 2AcA 



Pi 



C^{xt) - In 



ml 



+ A 



(19) 



The first term in this expression is the well known leading term, while two of the subleading 
terms are proportional to mf and carry the same weak phase as the leading term. In these 
contributions we have the two unknown matching functions C12 and C'12] however, the running 
gives the large logarithm which in this framework is assumed to dominate. Furthermore, the 
renormalization group reproduces the well known result for the term proportional to m^, which 
will modify the mixing phase. Note that this effective Hamiltonian is given in terms of local 
operators Pi. We shall estimate the matrix elements of this effective Hamiltonian using naive 
factorization at the scale /i ~ mc- 

From this effective Hamiltonian we can obtain the correction to the mixing phase. Numer- 
ically we find 



Aim 



M12 " 


1 ml 


\Vcd\ 


\Vct\ 




Coixt) ml 


\Vtd\ 


\Vtb\ 



\Vcd 


\Vcb\ 


\Vtd 


iVtbl 



sin(2/3) cos(2/3) + 2 C^ixt) - In 



(sin 13 — sin(2/9) cos [3 
= -(4.48 ±2.55) ■ 10-^ (20) 



At tree level, the next order in the 1/mb expansion consists of operators of dimension-9, 
which are six-quark operators (cf. e. g. jHj). They originate from diagrams like {ip = u,c) 




The contributions from these operators can be brought into the form 



(5°|C)6quarks|5°) ~ {B'\{Kd) {Kd){uu)\B'') 



m 



mt 



ml 



(21) 



where ^ is defined as 

{B'^\{hJ)(hJ)(uu)\W) 



1 + e (22) 



{B%Kd){Kd){cc)\BO) 

The GIM mechanism as well as the OPE guarantee, that ^ is of the same order as the m^/m^ 
terms. Altogether, the dim-9 operators are suppressed with respect to the dim-6 operators by 
a factor of AqQf^/ml, as could have been guessed from the OPE. Hence they are a negligible 
contribution, at most of the absolute order of 10~^. 

Another class of corrections are the 0{as) QCD corrections which are known for most of the 
processes at the two- loop level [12]. For the box diagrams they have been already calculated 
some time ago (see the references in However, one may use the effective field theory 

to resum large logarithms of the form ln(M^//i), where is a hadronic scale. Since there 
are large logarithms already at order (q!s)° the situation is similar to the one in the transition 
s —* dii, which has been discussed in [TT)] . 

The result given in (fT^ does not resum the logarithms of the form ln(Mvi///i). For the 
case of the AS* = 2 effective Hamiltonian, the next-to- leading result has been given in jllj . 
however, in our case the situation is slightly different due to the fact that the mass of the 
bottom quark sets a large scale and thus a matching to an effective theory with a static b quark 
at the scale fi ~ nib is possible, which is still perturbative. The running below the scale down 
to the scale rric resums logarithms of the form In m^/m^; however, these logarithms are not large 
and hence the resummation is not really needed. In principle one could run perturbatively even 
below the charm mass, yielding a result like ()19|) in terms of local operators, but the running 
below mb is a very small effect which we shall neglect in the following. 

Working at one loop requires to include the one-loop correction to the AB = 1 effective 
Hamiltonian as well as the mixing among the local AB = 2 operators of dimension 8. However, 
since we are only interested in the contributions that modify the relation between sin(2/5) 
and Sb^j/^Ks we shall simplify the discussion by neglecting the mixing among the AB = 2 
operators of dimension 8, since only the single dim-8 operator proportional to will contribute 
to this effect. 

The running of the AB = 1 effective Hamiltonian forces us to introduce operators Tj with dif- 
ferent color combinations. It is well known that for the AB = 1 operators the renormalization- 
group evolution is diagonalized by the combinations jTTj 

{{bLl.qLMl'dL))^ = ^ [{hLl,qi,L){q;,Ll%,L) ± {kLl,qj,L){q;,Ll''d,,L)] • (23) 
Thus it is convenient to introduce the non-local operators in the form 



( {bLl,qL){x) {q'^l^d,) {x))^{ ihr q'L) (0) (qLlM (0)) ^, 



(24) 



where a and a' may take the values +, — . Ignoring the mixing among the local dim-8 operators 
we choose as the basis 



'-^ IV 3 ) cc ' cc ' cc ' cn ' cn ' c« ^ ''- uu ■> uu ■> uu ) " 
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For the case at hand the relevant contribution is the mixing of the operators from the time- 
ordered products into the local operator Q^, for which we also keep the mixing with itself. 
The anomalous dimension matrix, including tree level and this restricted set of corrections, 
becomes 



7 = 
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IT 
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3 

167r2 
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2^ 
-^^ 

TT 

-A^ 












167r2 
1 

1fi7r2 


3 

327r2 
1 







2^ 

TT 

0-^0 
-4^ 

TT 







327r2 
1 

3?!,r2 

















2^ 

TT 

0-^0 
-4^ J 



(25) 



In terms of the "diagonalized" coefficients 



C'^'^ — C22 + Cii + C12 + C21 

C = C22 + C\\ — Ci2 — C21 

C+ = C22 ~ Cii + C12 — C21 

C + = C22 ~ Cii — C12 + C21 



(26) 



the matching conditions (i.e. the starting points for the renormalization-group evolution) are 



C++(M^) = A 



w) 



Ct-{M^) = 2Xl 
C.V(^^) = 2A^, 



\2 



c;-iM^) = Xi 

C~{M^) = 2XAu: (27) 



where we combined H — and — h. 

The sohition of the renormalization-group equation (at one loop) for the Wilson coefficients 
of the AB — 1 operators is well known {j^" — j'^^ ^^(/x)/?!, — 2, —1, —4) and yields 



C^5' if,) ^ C^Z' {Mw) 



as(/i) 



as{Mw)J ' ^^^^ 

where Po = 11 — 2iV//3 is as usual the coefficient of the one- loop QCD beta function. Due to 
the special structure of the anomalous dimension matrix, we obtain the renormalization-group 
equation for the coefficient C3 of the operator 



o^sil^)^ / ..X AcAt 



TT 



167r2 



Qs(/i) 
asiMw) 



-27V/30 



(29) 
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where r^"^'"* ' = 3,-2,1 originates from the first column of ()25j] . Note that there is a 
factor two for H — due to the equal contribution from H — and — h. Furthermore, the fiavor 
dependence of the time-ordered products' mixing into the quasi-local operators appeared within 
the structure + | ■ 2A„Ac = — AcA*. The factor of two here stems from the fact that the uc 
and cu fiavor combinations give the same contribution. 

The equation for C3 can be solved by standard methods; the solution of the homogeneous 
differential equation is 



W) 



-273 //3o 



By setting C^lfi) = c{fi)C^\fi) we get the solution 



(30) 



w 



-273//30 



C3{Mw,xt) + 



AcAt 
8^ 

1 



2(73-7?) 



/3o 



as(/i) 



as(M, 



w 



■ (31) 



We shall consider only the evolution from down to mi,, since at mj, we would need to 
consider again a different set of operators, including static quarks for the b and later also for 
the c and their renormalization ^Hj- While this can be done in principle, the corresponding 
logarithms ln(mf,/mc) and ag ln(mc/ /i) are smaller than the ones from the running from Mw 
to mft and we shall include here only the leading term. 

We obtain for the QCD corrections at /i ~ the explicit formula 



Aim 



M12 " 


1 ml 


\Vcd\ 


\Vcb\ 


IM12L 


Co{xt) m^ 


\Vtd\ 


\V4 



-273//30 



\Vcd\ 


\Vcb\ 


\Vtd\ 


\Vtb\ 



27r 

To 



2(73-7^) ' ' 



sin(2/?) cos(2/?) 



as(/i) \asiMw) 



/3o 



(sin/5 -sin(2/3) COS /5 ) ^, (32) 




which turns into the simple case (cf. (f^ ) as as{Mw) — > 0. Numerically, one obtains at = mft 

M12 

QCD 



Aim 



I Ml 



12 



-2.08-10- 



(33) 



which has to be compared to 



Aim 



Tm- 



12 



12 



-3.00 ■ 10" 



(34) 
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without QCD corrections. This indicates that QCD can reduce the absolute value of the 
correction to sin(2/?) coming from mixing by about 30%. Although we have not included the 
mixing among the dim-8 operators, we still take this as the size of the QCD corrections to be 
expected. We conclude that these contributions are safely below the percent level. 



4 Corrections to the decay 



The operators contributing to the decay to J/'^Ks are all of the flavor structure {bq){qs) 
where at leading order we have from the current-current operators q = c,u. From these two 
operators g = c is Cabibbo-favored over q = u. Furthermore, also the penguin contributions 
are dominated by the charm quark, such that all contributions to the decay carry the same 
weak phase and hence one expects a very small direct CP asymmetry Note that the top 
quark contribution has been integrated out already at the weak scale for both the QCD and 
the electroweak penguins. 

However, looking for small deviations we have to study the small terms carrying different 
weak phases which is the (bu){us) contribution. Aside from the QCD penguins also elec- 
troweak penguins become important once small corrections are studied ^H]- We identify the 
corresponding matrix elements of the two different contributions - tree and penguin - to the 
effective Hamiltonian in the numerator and denominator of the ratio of amplitudes 

I ^— i-y 

(35) 

where %s{b qqs) is the sum over the operators (multiplied by their Wilson Coefficients) with 
the quark content {b — >■ qqs). It is convenient to define the ratio 



A(SO - 






%s{b - 


^ccs)\B^) + {J/^Ks 


Zs{b - 


-^uus)\B'') lU^e-'-r 


A{B^ - 






Zs{b - 


ccs)\BO) + {J/'^Ks 


%s{b - 





r :- 



{J/^>Ks\%s{b^uus)\B^) 



{J/^Ks\Zs{b-^ccs)\B^) 
as an expansion parameter. In this way we get 



1 + re 



\Vub\ 


\Vus\ 


\Vcb\ 


\Vcs\ 



0.0203 ± 0.0066 



-17 



-17 



1 — 2 ir sin7. 



(36) 



(37) 



A{BO Jl-^Ks) 1 + re^ 

where r]j/^Ks = — 1 is the CP eigenvalue. 

The main obstacle to obtain a reliable quantitative estimate is the evaluation of these 
hadronic matrix elements. Some time ago the so-called BSS mechanism [201 has been sug- 
gested, where the ww-loops are evaluated perturbatively, assuming a sufficiently large momen- 
tum transfer through this loop. In fact, this approach has been supported recently by QCD 
factorization [21], indicating that the loop is indeed perturbative. 
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The evaluation of the loop requires to insert a typical momentum transfer k"^ passing through 
the up-quark loop [20] '■, furthermore, the loop also depends on the typical scale /i of the problem, 
which will be the mass of the h quark. Since the loop involves only scales well above the hadronic 
scale, one obtains again local operators and thus one may use as an estimate 



V2\ 



37r 



where the first term originates from the electroweak penguin, the second one from the QCD 
penguin. 

The remaining problem is to estimate the color-singlet and the color-octet matrix element. 
We shall take a simple-minded approach and estimate the sizes of these matrix elements from 
the total rate of the decay J /'^Ks- From the measured life time, r = (1.537±0.015) ps 

[22], and the branching ratio BR(5° J/'^Ks) = (4.25 ± 0.25) ■ 10"^ ^ one can calculate 
the matrix element for the cc-contribution by taking the square root 

I {J/-^Ks\{sb){cc)\B^) lexp. = (8.36 ± 0.66) • 10^ MeVl (39) 

The matrix element with the charm pair being in the octet is estimated by splitting the 
effective Hamiltonian in singlet and octet contributions. In naive factorization only the singlet 
piece survives which, however, yields a rate roughly a factor three to four too low compared 
to experiment. The difference comes from nonfactorizable contributions which we shall ascribe 
completely to the octet term. To obtain an estimate of this matrix element from the measured 
total rate we also need the relative phase of the two contributions, which we take from QCD 
light-cone sum rule estimates in [23 (see also [21] )• The relative phase turns out to be small 
and so we may simply add the two parts. 

For the singlet matrix element calculated by means of factorization we get 

I {J/^Ks\{sb){cc)\B'^) Ifact. = (3.96 ± 0.36) • 10^ MeVl (40) 
Inserting for the Wilson coefficients the values^ 

= Ci + -^C2 = 0.10 ±0.03, C^^) = 2C2 = 2.24 ± 0.04, (41) 

we get as an estimate for the octet matrix element 

I {J/^Ks\{sT%){cT''c)\^) I = (1.97 ± 0.64) ■ 10^ MeVl (42) 

^We take the Wilson coefficients at Leading-Log order at the scale mi,; unfortunately, the scale dependence 
is still large at leading order, but we are aiming only at an estimate. 
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Taking the scale to be mt and momentum transfer as mj/^ gives (the two contributions are 
the electroweak and the QCD penguins, respectively) 

r = [-(0.16 + 1.27) ±0.66] ■ 10"^ • - In (^^^ + ivr^ 

Re [r] = (-3.62 ± 1.55) ■ 10-^ Im [r] = (-4.48 ± 1.92 )■ 10"^ (43) 

We have to point out that the estimates of these matrix elements are extremely difficult and 
hence quite uncertain; the electroweak penguins have been estimated in a recent paper with a 
similar approach 



5 Effects on the CP asymmetry 

Now we are ready to collect all elements for the CP violation terms within the asymmetry. 
The quantity |A|^ is close to unity, the deviation being a small quantity of the same order as r. 
Thus we write |A|^ =: 1 + 6 and get 



C 



2' 



Therein 6 is given by 



6 = AA + e = 4 Im [r] sin 7 — Im 



12 



M- 



12 



(44) 
(45) 



where A A is the deviation of the [A/^l^ from 1 and e has been defined in (0]). 

For the ratio between the off-diagonal width and mass matrix elements one can take e.g. 
the calculation in |26p: 



Tm 



12 



Ml 



12 



An mi 

Im 



(5.18 + 2.96) ■ 10' 



from which we obtain a numerical value for 5 

5 = -(1.02 + 0.75) • 10"^ 



(46) 



(47) 



Including all small corrections, the time- dependent CP asymmetry in J /'^Ks has to 

be fitted to 



acp(^) = — (sin(2/3) + AS^-^j/^i^g) ■ sin(Amt) — - cos(Amt) 

-\ sin^(Amt) sin^(2/5) + sin(4/5)— — sin(Amt). 



(48) 



Besides the correction to the leading sin(2/5) term and the cosine term there is now a small 
constant contribution to the CP asymmetry as well as terms proportional to s\.T?{Amt) and 
t sin(Amt). 

•^Here we are assuming the solution for (3 in accord with the unitarity triangle. 
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The correction to the mixing-induced CP violation is (note that the term containing the 
imaginary part of r cancels in the expression for the correction) 



ASb^j/^ Ks = 2Im[r] sin 7 sin(2/5) — Aim 
= 2 sin 7 Re [r] cos(2/?) — Aim 



Ml 



12 



IM12I 
M12 



+ 2 sin 7 Re 



M- 



12 



I Ml 



12 



M 



121 



(49) 



With the most recent value of sin(2/9) = 0.736 ± 0.49 P, the corrections in ASb-^j/^ii Ks 
have the size of 



Aim 



M 



12 



IM 



12 



(2.08 ± 1.23) ■ 10" 



and sum up to 



2sin7Re[r]cos(2/3) = -(4.24 ± 1.94) ■ 10" 



ASB^j/'i,Ks = -(2-16 ± 2.23) ■ 10" 



(50a) 
(50b) 

(51) 



which is a correction of roughly a third of a per mil with respect to the measured value for 
sin(2/3), but due to the large uncertainty could also be much smaller. 

From the mass difference as measured quantity and the approximate calculation for 



Ar 



3tt ml 



1 - 



8 ml 




\Vcd\ 


3 ml 


Vtb\ 


\Vtd\ 



COS j3 



AMf 



(52) 



one can determine the term leading to a linear dependence on t in 

Ar ^ - (1.773 ± 0.249) • 10"^^ MeV 
^ - (2.694 ± 0.378) ■ 10"^ ps"\ 

Hence, this term has a typical size of the order of 

sin(4/3) ^^7^° ^ -(1.03 ± 0.15) ■ 10'^ , 



(53) 



(54) 



again a very small contribution, which, in addition, will not show up at small times since its 
time dependence is (t/r^o) sin(Ami). 



6 Discussion 

We have reinvestigated the well known fact, that the mixing-induced CP asymmetry of 
J/^>Ks provides us with a very clean measurement of sin (2/5). Already in the original paper 
[3j it was argued that in the standard model the contamination from "wrong" weak phases is 
tiny in this decay, but in the meantime the measurements became so precise that we considered 
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it worthwhile to attempt again to quantitatively analyze these small standard-model contribu- 
tions. In particular, we have used an effective field theory ansatz for the analysis of AB = 2, 
which has not been employed before to this process. 

Our motivation was twofold. First of all, it has been argued that the CP asymmetry could 
give a hint to new physics, which has been analyzed generically in 0. It turns out that the 
general picture conjectured in |lj is supported by the present analysis. Secondly, the B factories 
are doing very well and produce a large amount of data. Already now the measurement of the 
CP asymmetry in the gold-plated mode is a precision measurement with uncertainties at the 
level of less than ten percent. In the near future this measurement will improve further. 

We have shown that the corrections to be expected in the standard model can partially be 
calculated systematically, namely the part originating from corrections to the mixing. Unfor- 
tunately, the second contribution, which is the one from the decay matrix element, is much 
harder to access, so we still cannot obtain a reliable estimate for the corrections to the CP 
asymmetry of the gold-plated mode. This situation could be improved once data on the decay 
Bs J/'^Ks becomes available, since some of the uncertainties could be eliminated using 
these data |2ZI- However, currently one has to use the methods that have been proposed by 
different authors and one can infer that the corrections will be very small, in the range of a few 
per mil. At least we can conclude from our analysis that there is still room for new physics in 
this observable, since a deviation from the standard-model prediction at the level of percents 
(which is not yet the experimental accuracy) would indicate the presence of new physics. 
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